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Abstract 

We present a consistent calculation of half-off-shell form factors in the pion- 
nucleon vertex and the nucleon self-energy. Numerical results are presented. 
Near the on-shell point the pion-nucleon vertex is dominated by the pseu- 
dovector coupling, while at large nucleon invariant masses we find a sizable 
pseudoscalar admixture. 



Key Words: Few-body systems, Off-shell form factors, Nucleon self-energy, Pion-nucleon 
vertex. 

1996 PACS: 13.75. Gx, 14.20 Dh, 21.45. +v 

February 9, 2008 



1 



Typeset using REYTf^X 



I. INTRODUCTION 



The structure of hadronic vertices, usually parametrized in terms of form factors, is im- 
portant in much of nuclear physics. The form factors may depend on the different invariants 
that can be constructed. In nucleon-meson or nucleon-photon vertices one often considers 
only the dependence on the momentum-squared of the meson or photon. In the present 
paper we will consider so-called off-shell form factors where the dependence is studied on 
the momentum-squared of one of the nucleons involved. 

Off-shell form factors are an ingredient in the description of physical processes. For 
example, nucleon-photon off-shell form factors have been shown to be important in models 
for proton-proton bremsstrahlung and virtual Compton scattering irNN and other 
nucleon-meson form factors with an explicit dependence on the momentum of one or both 
nucleons have been used in models for NN || and ttN scattering, pion photoproduction 
H and vector meson production in nucleon-nucleon collisions [TO]. In these models, the 



form factors have been phenomenologically parametrized, with the parameters adjusted to 
fit experimental data. 

The off-shell structure of the nucleon-photon vertex [O-O], and the nucleon-pion vertex 



lll)14|15| has been studied before. In particular, dispersion relation techniques are used in 



Refs. [n,r2,H], whereas the models of Refs. |H| are based on a perturbative dressing of the 
vertex with one-meson loops. In this work we investigate the pion-nucleon coupling in a field- 
theoretical model which is inherently non-perturbative and is based on the Schwinger- Dyson 
equation, considering loops to all orders. The nucleon self-energy and the pion-nucleon vertex 
function are both calculated in a consistent framework. 

In general, off-shell form factors and the functions parametrizing the self-energy are 
complex functions, where the imaginary parts are related to open multi-particle chanels of 
which the pion-nucleon channel will be the most important. Our approach is based on the 
analytic structure fl6l , |Tll of the nucleon self-energy and the off-shell ttNN vertex, which is 
exploited by the use of dispersion relations. The imaginary parts of the form factors and 



the self-energy are calculated from Cutkosky rules [[L7]]. To make this procedure tractable, 
we consistently neglect contributions to the imaginary parts from the multi-pion thresholds. 

For the course of this paper, we are interested in vertices with one off-shell nucleon, which 
contain two independent form factors. The convergence of the loop corrections is insured 
through the introduction of a "cut-off" function (the initial form factor). We obtained the 
interesting result that the widths of the converged form factors have an upper bound. 

Solutions of the Schwinger-Dyson equation have been presented in the past (see, e.g., 
T8| , |T9ll and a recent paper pPf). There, a usually adopted approximation consists in assuming 



the same spin structure for the dressed and the bare vertex. The dressing of the vertex is 
thus parametrized in terms of a single form factor. In the present work we have released 
this condition and found a strong dependence of the spin structure of the vertex on the 
off-shellness involved. 

Form factors are usually interpreted as representing the features that are not included ex- 
plicitely in a particular model for a physical process, and as such, should be built consistently 
with the kind of models in which they are intended to be used. The form factors considered 
in the present paper are primarily designed for usage in a K-matrix model for pion-nucleon 



scattering, pion photoproduction and Compton scattering off the nucleon [21,22]. Since in 
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such a model the one-pion production channel is included explicitely, only real form factor 
should be used there (at least below the two-pion threshold). This aspect is elaborated on 
in Section 1111 C. 



In any model, treatment of off-shell three-point vertices should be linked with treatment 
of higher-point vertices, because a redefinition of the nucleon field can change off-shell de- 
pendence of the former in favour of presence of the latter. The observables are oblivious to 
the representation of fields (this result is known as the equivalence theorem) examples 
of which can be found, e.g., in Refs. p4]- p6[| . In the present model, higher-point vertices are 
excluded at all stages of the calculations, and the discussion is carried out solely in terms of 
off-shell form factors in the tcNN vertex. 

The paper is organized as follows. In Section II the general structure of the off-shell 
ttNN vertex is discussed. Our model is described in detail in Section III. At present we 
limit ourselves to the inclusion of one-pion-nucleon loops only, for which numerical results 
are presented in Section IV. 



II. STRUCTURE OF THE irNN VERTEX 

The tiNN vertex operator is the sum of all connected Feynman diagrams with one 
incoming nucleon (carrying the momentum p), one outgoing nucleon (p') and one pion 
(q = p — p'), with the propagators for the external legs stripped away. The most general 



form compatible with Lorentz covariance and isospin invariance reads |£7| 



r a ( P ',P,q) = t^g^MW^g^M + *—^G,{p'\p\?) 

m m 



p" -m A-m /2 2 2 



^—^G 4 (p> 2 ,p 2 ,q 2 )), (1) 
mm ' 

where m denotes the nucleon mass and r a ,a = 1,2,3, are the isospin Pauli matrices. The 
form factors Gi depend on the three Lorentz scalars, p' 2 ,p 2 and q 2 . Usually the situation 
is considered in which both nucleons are on the mass shell, i.e. p' 2 = p 2 = m 2 , and only 
Gi(m 2 , m 2 , q 2 ) enters in Eq. ([!]). In this paper we consider a different situation in which the 
pion and only one of the nucleons is on the respective mass shell, p' 2 = m 2 and q 2 = fi 2 , 
where /i denotes the pion mass. Such a vertex is conventionally called the half-off-shell 7iNN 
vertex, and it contains so-called half-off-shell form factors. 

If the operator of Eq. ([!]) works on the positive energy spinor u(p') to the left, the last two 
terms in Eq. ([TJ) vanish due to the Dirac equation, u(p') ft = u{p') m, and the vertex contains 
only the form factors Gi(m 2 ,p 2 , /i 2 ) and G2(m 2 ,p 2 , fi 2 ). Similarly, if the initial nucleon is on- 
shell, only the form factors Gi(p' 2 , m 2 , /i 2 ) and G^{p' 2 , m 2 , /i 2 ) are left. Charge-conjugation, 
space-inversion and time-reversal symmetries allow to relate these form factors, 

G 1 (p 2 ,m 2 ,n 2 )=G 1 (m 2 ,p 2 , f i 2 ), G 3 (p 2 , m 2 , ^) = G 2 (m 2 ,p 2 , /i 2 ). (2) 

Hence, one can consider only the vertex with the outgoing on-shell nucleon. Omitting the 
trivial arguments in Gi, 

T a (m,p,/j) = T a T(m,p,jj) = T a ^[Gi(p 2 ) H G 2 (p 2 )), (3) 
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where the notation T a (m,p, p) implies that q 2 = p 2 and u(p') p" = u(p') m in all expressions 
for this vertex. Along with Eq. (|3|), we will use another form for the half-off- shell vertex, 

T a (m,p,p) = r al \G PS {p 2 ) + ^^G PV (p 2 )), (4) 

where 

Gps{p 2 ) = G 1 (p 2 ) - 2G 2 (p 2 ), G PV {p 2 ) = 2G 2 {p 2 ) , (5) 
denote the form factors corresponding to the usual pseudoscalar and pseudovector couplings. 

III. DESCRIPTION OF THE MODEL 

The model for the form factors is based on a non-perturbative dressing of the vertex with 
pion loops as represented graphically in Fig. 1. The nucleon self-energy is calculated 
self-consistently using the Schwinger-Dyson equation E3] with the dressed vertex. This 
can be expressed in terms of a system of integral equations for the dressed wNN vertex 
T a (p',p,q) and the dressed nucleon propagator S(p), 

f d 4 k i 

r a (m,p,//) = r° a (m,p,n) -ij j^-^(Tp(m,p' + k,k)S(p' + k) 

xT a (p' + k,p + k,p) S(p + k) T p {p + k, p, -k) D{k 2 )) , 
S(p) = S°(p) + S(p)J:(p)S (p), (6) 

S(p) = -ij ^{Y^p + k^Sip + ^Vlip + k^-^Dik 2 )) 

-(Z 2 - l){j>-m) - Z 2 5m, 

where D is the pion propagator, S° the free propagator of the nucleon and r°(m,p, fi) the 
bare ttNN vertex. The last two terms in the equation for the self-energy are part of the 
renormalization procedure and will be discussed later. The dressing is non-perturbative 
since the dressed vertex and propagator appear also on the right-hand side of the equations. 
As is well known, such a procedure suffers from divergences. In addition, we are interested to 
build a model for half-off- shell vertices while the right-hand sides of Eqs. @ include vertices 
with all external legs off shell. To circumvent these two problems we have applied a solution 
procedure based on the use of dispersion relations and a regularization method as outlined 
in the following sections. 

A. Solution procedure 

Hereafter we shall denote half-off-shell vertices as r a (p), dropping the trivial parameters 
for brevity. As noted in the above, Eqs. (0) require the knowledge of the full off-shell vertex. 
Using the analyticity of the form factors and the self-energy this problem can be bypassed. 
The imaginary parts of the form factors can be obtained by applying Cutkosky rules to the 
integrals in Eqs. (|6]), and we reconstruct the real parts through the application of dispersion 
relations implemented in an iterative procedure. 
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As stated in the introduction, we have excluded four- or higher-point vertices from our 
model. Furthermore, in applying Cutkosky rules we shall only include the channel with 
the lowest threshold, i.e. the one-pion nucleon channel. This results in the contributions 
depicted in Fig. 2. 

The solution procedure can now be explained best by going in some detail through one 
complete step of the iterative procedure. From the n th iteration we have obtained the form 
factors Gi 2 {p 2 ) which define the vertex through Eq. (^) and the self-energy functions A n (p 2 ) 
and B n {jp 2 ) which define the dressed propagator, 

(S n (p)y l = Z 2 ty- m B ) - [A n {p 2 )i> + B n (p 2 )m] , (7) 

where the term in square brackets is the loop contribution S^(p) to the self-energy. The 
parameters Z 2 and mg = m — 5m are renormalization constants as defined in Sec. (|IHB|) . 



The imaginary parts of the form factors and the self-energy functions arise from the 



pinching-pole term in the loop integrals which can be evaluated using Cutkosky rules |T7 
This contribution is labelled by the subscript I. For the self-energy one has: 

= -r:>) r^>) , (8) 

where the subscript R denotes the vertex calculated using only the real parts of the form 
factors Gi2- The Dirac conjugated vertex is denoted as T. The remaining integral can be 
written as 



'■pole 



The explicit form for I po i e {p) is given in Appendix A. The imaginary parts of 
A n+1 (p 2 ) , B n+1 (p 2 ) can now readily be written using Eq. (J7J). 



The real parts of the self-energy functions are calculated via the dispersion relations [16 



7T Jw 2 th p' 2 — p 2 

and similar for ReB(p 2 ). Here w t h = (m + [/,) is the one-pion threshold, and V denotes the 
principal value integral. 

The pole term in the loop integral for the vertex reads: 

Kj 1 (p) = J P oMK,r(p), (ii) 

with 

J P oUp) = - ^ I d A k T n R (p' + k) S n R + \p' + k) T n R (p' + k){p ! + k ! + m) 

x8{(p + k) 2 - m 2 ) 6(p + ho) 5(k 2 - fi 2 ) 9(-k ) , (12) 



where the integral Eq. ([12]) is independent of p', the momentum of the outgoing on-shell 
nucleon, as shown in Appendix B. The factor (—2) in Eq. (12) comes from commuting the 
isospin matrices. The propagator S^ +1 in Eq. (|T2"D contains only the real parts of A n+1 and 
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B n+l which have been calculated by virtue of Eqs. (^|,|9|JT0D . Casting r""j (p) in the form of 
Eq. ([3]), the imaginary parts of the form factors Gi^ 1 (p 2 ) are found. 

To construct the real parts of the form factors we take advantage of their analytical 
properties [11], 



Re GTV) = C<V) + V - fV . (13) 

7T Jwj h p -p 

The first term on the right-hand side of Eq. ( |T3"D derives from the equivalent term in Eqs. (^|). 
We use unsubtracted dispersion relations since convergence of the integrals can be guaranteed 
by the cut-off function introduced in Eq. (|18|). 

There are a few points that need special stressing here. 

• In calculating the imaginary parts for the (n+l) st iteration, we retain only the real 
parts of the form factors and self-energy functions from the n th iteration, to be consis- 
tent with the use in a K- matrix formalism as explained in Sec. ( [Ill C| ) . 

• Except the depicted cuts, any other kinematically allowed cuts of the loop diagrams 
(cutting through the blobs in Fig. (|2|)) would correspond to either picking up contri- 
butions of higher thresholds or considering four-point vertices, both of which would be 
inconsistent with the adopted solution scheme. 

• In cutting the self-energy loop diagram, dressed vertices at both sides of the cut prop- 
agator are taken into account. Seemingly, this is in conflict with the Schwinger- Dyson 
equation (see Eqs. (^) and Fig. 1), where the second blob would lead to double count- 
ing. However, the presence of the two blobs in the cut diagram is necessary to sum up 
all contributions from one-pion-nucleon cuts. 

• The present method of solution allows one to avoid dealing with the full off-shell 
vertices present Eqs. (|6|). Indeed, as can be seen from Fig. 2, we need only half-off- 
shell vertices throughout the iteration process. 

• To calculate the pole contributions of the loop integrals we have applied Cutkosky 
rules i.e. put the nucleon and pion lines in the loop integrals on their respective mass 
shells, as shown in Fig. 2. In the cut propagators, therefore, only physical masses 
appear. In particular, this implies that the dressing of the pion propagator does not 
have to be considered in the present approach. 

We take r°(m,p, /i) as zeroth iteration for the vertex (its precise form is specified in 
Eq. (|l"8D). At each iteration step, we utilize the dispersion relations Eqs. (|T3|JT0|) , where 
Im Gi i2 (p 2 ), Im A{p 2 ) and ImB{p 2 ) are calculated using Cutkosky rules. In this connection, 
the following remarks are in order. 

Analyticity. 

In principle, the use of the dispersion relations should guarantee that the form factors and 
the self-energy functions A and B be analytic in the complex plane cut from w 2 h to oo along 
the real axis. However, the actual imagianary parts calculated in the model contain also 
"unphysical" singularities of the function G°(p 2 ), see Eq. ([18|), regularizing the dispersion 
integrals. This, strictly speaking, invalidates the derivation of the dispersion relations. This 
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problem will be encountered for any non-constant function G°(p 2 ) unless its singularities are 
located along the cut. 

Unitarity 

In applying Cutkosky rules, an important qualification is that we neglect the contributions 
to the imaginary parts that come from the intermediate states including one nucleon and 
more than one pion. In order that unitarity should hold exactly, the imaginary parts must 
contain the contributions from all multi-pion thresholds. In the context of the present work, 
a rigorous account of, e.g., the two-pion threshold would require computing the imaginary 
parts of two-loop self-energy and vertex diagrams. Analyses of massive two-loop Feynman 
diagrams have appeared in the literature recently, including the dispersion relation approach 
(see, e.g., [29| and references therein). However, such calculations are rather complicated, 
and we found their application in the present model not feasible. 

One may argue that the real function G°, which is used presently, could be expressed 
in terms of a dispersion integral over a function, say F°, which would correspond to the 
discontinuities of the form factors due to all channels opening at higher thresholds not 
considered explicitely. If this were the case, the presently adopted procedure would be 
equivalent to adding an extra contribution F° to the imaginary parts derived from the loop 
diagrams. This could possibly account for the two problems just mentioned, but we did not 
pursue this direction. 



B. Renormalization and regularization 

The renormalized nucleon self-energy in Eqs. (^) or Eq. (|7]) can be written as 

S(p) = E L (p) - (Z 2 - 1)$ - m) - Z 2 5m. (14) 

The first term in Eq. ( |TS ) is the contribution of pion loops while the last two terms come 
from the counterterms in the Lagrangian as part of the renormalization procedure. 

The construction of the counterterms is based on the usual renormalization procedure 
30j as explained by the following example. The Lagrangian, written in terms of the "bare" 



fields, masses and coupling constant, is 

£ = hd v <f> B d»(f> B - /i| 0|) + $ B (i p - m B )ll>B ~ 7T~ i^Bl 5 (P<f>B)i>B. (15) 

Defining the renormalized nucleon field ift = Z^ 1 ^ 2 ipB> the renormalized nucleon mass m = 
tub + 5m and the constant //(2m) = gBZ 2 /(2m B ), Eq. (|l~5|) can be reformulated as 

C = ~ (d u <p d u 4> - fi 2 ^)+^(ip-m)ij 

~ (if- V¥ (s^>)V> - Z 2 5m^ - (Z 2 - 1) ^(i ft - mjip) . (16) 

Because we encounter only cut pion lines during the iteration procedure (see Fig. 2), we 
need only the pole contribution of the pion propagator. Since the pole properties of the 
renormalized dressed propagator coincide with those of the free one, the pion field and mass 
need not be renormalized in our approach, = jj? = n B . From Eq. (0) and Eq. (|j~4j) it 
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can be seen that the renormalization constants Z 2 and 8m can also be interpreted as real 
constants which can always be added when the real part of a function is determined from 
the imaginary part via a dispersion relation Eq. (|T0|). 

The coupling strength / and the renormalization constants Z 2 , 8m are determined by 
fixing T a (m, m, /i) and the pole structure of the propagator S(p), 

5 _1 (m) = 0, 

Res[S(p),j = m] = l, (17) 
u{p') r a (m,m,ii) u(p) = u(p') r„7 5 gu(p), 

where the last equation can be reduced to Gi{m 2 ) = g, the physical pion-nucleon coupling 
constant (we take g = 13.02 ||31|| ). The left-hand side of this condition is calculated at the 
kinematically forbidden point, where all the external legs of the vertex are on-shell. However, 
this is of no harm for the renormalization prescription. We could choose any convenient 
renormalization point as long as the form factors calculated at that point are real (see, e.g. 
PDfl , where the freedom of the choice of a renormalization procedure is discussed in general) . 

In the context of the iterative procedure described in the previous section, the constants 
Z2 and 8m are chosen to provide the correct pole properties of the converged propagator. 
This implies that the pole location and residue of the propagator are off in the course of the 
first few iterations. To check that this feature is immaterial for the final result, we applied 
also another solution procedure. Its main difference from the one outlined above is that 
the renormalization of the propagator is done at each iteration step, insuring the correct 
pole properties at any iteration. We found that both methods lead to identical results for 
the converged vertex and propagator. The reason for this is that when convergence has 
been reached, a non-perturbative solution of Eqs. (^) (under the provisions which have been 
discussed) has been obtained. The intermediate steps in the iteration procedure at this point 
are an uninteresting technical detail. 



The loop integrals, or rather the dispersion integrals Eqs. ( |I3| , |10|) , will diverge unless 
a regularization is applied. As part of the regularization procedure, we introduced a form 
factor G°{p 2 ) (also called the cut-off function) for the bare nNN vertex, 

Tl(m,p^)=r a ^-±^G°( P 2 ), (18) 

in terms of Eq. (^). The function G°(p 2 ) = G° PV (p 2 ) is normalized to / at p 2 = m 2 and must 
fall off sufficiently fast at infinity to provide convergence of the integrals. In the numerical 
example discussed later we used two different functions G°(p 2 ), see Eqs. 



The cut-off function is a phenomenological input of the model. A self-consistent proce- 
dure to construct meson- nucleon form factors was presented in Ref . |32| , where both nucleons 
in the vertex are on-shell and the meson is off-shell. There, no phenomenological form factor 
was needed. We were not successful in implementing a similar approach for half-off-shell 
form factors in the pion-nucleon vertex, nontrivial solutions of the relevant equations for the 
self-energy and the 7rNN vertex did not seem to exist. 
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C. Consistency with a K-matrix approach 



One motivation for the present model is the construction of form factors and self-energies 



which can be applied in a K-matrix approach to ttN scattering [EUpIJ. We outline the K 



matrix method (details can be found in plf ) and in particular address the double-counting 



issue: by considering only the real part of the form factors and the self-energy functions 
on the r.h.s. of Eqs. (|3| |lT| , |T2|) , we avoid double counting when the calcualted vertex and 



propagator are used in the K-matrix approach. It should be emphasized that only the 
one-pion threshold discontinuities are taken into account in both the K-matrix approach in 
question |^l],^2j and the present model. 



The Bethe-Salpeter equation for the irN scattering amplitude T can be written in the 
operator form, 

t = v + vgr. (19) 

Here, V is the sum of all irreducible diagrams describing the scattering, and Q is the free 
ttN propagator. Q can be decomposed as the sum of the on-shell contribution i5 which is 
imaginary (according to Cutkosky rules), and the off-shell part Q v which is real, 

g = g p + is, (20) 

where 5 implies that the corresponding intermediate nucleon and pion are taken on their 
respective mass shells. The K-operator is introduced by the equation 

K = V + VG P K. (21) 

Combining the last three equations yields the T-matrix expressed in terms of the K-matrix, 

T = K + Ki5T. (22) 

This is the central equation used in a K-matrix approach and can schematically be written 
as T = Kj (1 — iK). If K is hermitian, the scattering operator 5=1-1- 2iT will be unitary. 

We consider a simplified version of the K-matrix approach containing only nucleons and 
pions, with the kernel V chosen as the sum of the s- and u-channel tree diagrams. One way to 
construct the K-operator is to set K = V p2| , thereby assuming g p = 0, see Eq. (^I|). Then, 
by Eq. (^), the T-matrix will contain the loop diagrams in which only the cut nucleon and 
pion propagators will enter. Only by using dressed vertices and propagators one may take 
the K-matrix equal to the sum of skeleton diagrams solely. As implied by Eq. (pip, the form 
factors in these dressed vertices take into account real contributions due to the principal 
value g p . These are the real parts of the form factors discussed in the previous sections. If 
we kept both the real and imaginary parts of the form factors and the self-energy functions 
in Eqs. (j^ |Tl~l , |l2"l) , it would be inconsistent with the K-matrix approach. In particlular, the 



on-shell contributions iS would be taken into account twice for every nN propagator g. An 
exception are some one-particle irreducible diagrams contributing to the T-matrix. 

It is known that the nucleon and pion degrees of freedom are not enough for a realistic 
description of 7rN scattering (for example, the role of both the Delta resonance and the p 
meson is indispensable) [p|-|9l,pT|,[Z2]] . Since in our model we confine ourselves to the pion 
and the nucleon, no calculations for the nN scattering observables will be presented in this 
paper. 
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IV. NUMERICAL RESULTS 



Two sets of calculations were done, corresponding to the two following cut-off function 
GV), Eq. ©: 



GW) = f 

and 



(A 2 — m 



2\2 



(A 2 — m 2 ) 2 + {p 2 — m 2 ) 2 



(23) 



GW) = f e~ JL ^~ . (24) 

The functional dependence of G%p 2 ) is taken from Ref. ||, where it was used as an off-shell 
form factor in the nNN vertex. We define a parameter A 2 = p\ 5 — m 2 , where p^ 5 is the point 
at which G°(p 2 ) reduces by factor two comparing to its maximum value / (here Pq 5 > m 2 ). 
Then, for the calculations with the functions Eq. (p3[) and Eq. A 2 equals 



A?(A) 

and 



(P ™1 )2 - (\2 - m 2)2 (25) 

\ Vol 



Ajj(d) = In 0.5, (26) 

respectively. We find that the iteration procedure described above converges only if A < 
A c m 1.7GeV for G°j(p 2 ), and if d < d c ps 1.65 for G°/(p 2 )- The corresponding "critical" 
values for the half- widths can be inferred from Eq. (E3) and Eq. (p6|): A 2 (A C ) = 1.28GeV 2 
and Af J (d c ) = 1.33GeV 2 . Results of calculations are presented below for the following two 

cases: 

Case (I) Calculations with the cut-off function Eq. (p3[), where A = A c = 1.7GeV; 
Case (II) Calculations with the cut-off function Eq. (EMI), where d = d c = 1.65. 



As stated above, the constants /, and 8m are chosen to satisfy Eqs. (|17|) . The values of 
these constants for cases (I) and (II) are given in Table 1. 

The convergence was considered achieved at iteration m if all the results of iterations 
m+l,...,m-|-20 were identical to those of iteration m up to six significant digits. We could 
put a very strong convergence criterion since the computer program uses little CPU-time. 
With this criterion, convergence was reached after about 100 iterations. We mention that, 
for example, the self-energy after 10 iterations differs still quite noticably from the converged 
result. 

A comparison of results obtained with the two different cut-off functions show how these 
reflect in the final results. The non-perturbative aspects are stressed by comparing the 
results of the first iteration (basically a one-loop calculation) with those of the converged 
calculation. 
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A. Results for the half-off-shell form factors 



The imaginary parts of the form factors Gpy(p 2 ) and Gps(p 2 ) are shown in Fig. 3. The 
results of calculations for the two cut-off functions introduced in Eq. (|23|) and Eq. (|24|) are 
shown next to each other. For case (II) the tails of the form factors at large off-shellness are 
suppressed due to the exponential in the cut-off function. Independent of the choice of the 
cut-off function there is a marked difference in the results of the first iteration (dotted curve) 
and the converged results for the pseudovector form factor. The reason for this difference 
is the (small) pseudoscalar component of the final form factor. The converged and first 
iteration results for the pseudoscalar form factor differ much less, as can be seen from the 
bottom panels of Fig. 3. 

The real parts of the form factors are shown in Fig. 4. The top pannels show the 
pseudovector form factor Gpy{p 2 ) (the solid line) together with the zeroth iteration form 
factor G PV {p 2 ) (the dotted line) which equals the cut-off functions Eq. ( P3"D (left) and Eq. ( |2"4] ) 
(right). It is seen that the bulk of Gpy{p 2 ) is contained already in G PV (p 2 ), and only a 
small part comes from the loop corrections. This manifests itself also in the small difference 
between the constant / and the physical coupling constant g, as can be read from Table 
1. We conclude that, in the present model, the shape of the converged form factor Gpy 
depends strongly on the phenomenologically introduced cut-off function. The middle panels 
of Fig 4 give more insight in the role of the pion dressing. There, the real part of the 
pseudoscalar form factor Gps{p 2 ) of the first iteration (the dashed line) is shown together 
with the converged result (the solid line). Since the zeroth iteration vertex is chosen purely 
pseudovector (Eq. ([TJD), Gpsip 2 ) appears solely due to the dressing. Also shown is the 
difference Re(Gpy{p 2 ) — G PV {p 2 )) which is the dressing contribution to the real part of the 
pseudovector form factor (the dash-dotted and dotted lines for the first iteration and the 
converged result, respectively). Note that the deviation of the non-perturbative result form 
that of the first iteration is considerable for this quantity. The ratio of the real parts of 
the Gps{p 2 )- and Gpv{p 2 )- form factors is given in the bottom panel of Fig 4. It is small 
below the pion threshold (e.g., Gps(w 2 h )/ g is about 2.1 % for both cases (I) and (II)), but 
becomes larger at higher p 2 . Note that Gpy{p 2 ) decreases for case (II) faster than for case 
(I), whereas the behaviour of Gps{p 2 ) for the two cases is comparable. This explains why 
the absolute value of (ReG psip 2 )) / (ReG pvip 2 )) grows faster for case (II) than for case (I). 

We remark that admixtures of the pseudovector and pseudoscalar pion-nucleon couplings 
have been studied in the past in connection with the NN and ttN scattering processes, where 
the vertex has been determined by adjusting phenomenological parameters to fit data (see 
discussions in |]5|-|7|,|2T|) . In those calculations the admixture is assumed to be constant. 



Instead, the present results indicate that the ratio is strongly dependent on the momentum 
of the off-shell nucleon. Evidence for large pseudoscalar admixtures for far off-shell momenta 
has also been observed in calculations of pion-photoproduction 



B. Results for the self-energy 

The imaginary and real parts of the functions A(P 2 ) and B(P 2 ) are shown in Figs. 5 and 
6, respectively. The solid (dotted) lines are the converged (first iteration) results. One can 
see that these functions approach zero faster for case (II) than for case (I). Of course, this is 
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entailed by the softer behaviour of Gp V (p 2 ) for case (II) as opposed to case (I) (see Fig. 4). 
The difference between the converged results and those of the first iteration is substantial, 
especially for the function B(p 2 ). The on-shell value E^m) = m(A(m 2 ) + B(m 2 )) = Z 2 <5m is 
negative and equals —48.8 MeV for case (I) and —47.2 MeV for case (II), see Table 1. Please 
note that if we had chosen smaller values for the cut-off, these self-energy corrections would 
have been less. For comparison, we mention that the contribution to the nucleon mass shift 
from one-pion loop calculated in baryon chiral perturbation theory yields a nucleon mass 
shift of about -15 MeV 0. 

Having obtained A(p 2 ) and B(p 2 ), one can find the spectral function of the self-energy 
T(uj) from Eq. (|A9|). Fig. 7 shows the spectral function for the two cases (the upper and 
lower panels). The dotted and dash-dotted lines are, respectively, the first iteration and the 
converged spectral function found in our model. In spite of the fact that ImB 1 (p 2 ) differs 
considerably from ImB(p 2 ), having the opposite sign at some momenta-squared (see Fig. 5), 
the spectral function remains positive for all iterations (as it should). 



The spectral function of the nucleon propagator was recently considered in Ref. f20 
whose approach is, however, different from the present work. In particular, there the vertex 
was not calculated consistently with the nucleon propagator. 



V. CONCLUSIONS 

We have presented a solution procedure for the Schwinger-Dyson equation to obtain 
consistently the nucleon self-energy and the half-off- shell pion-nucleon vertex. Retaining 
the non-perturbative aspects of the equations is important. We observe a large difference 
between the simple one-loop results and those of the converged procedure. As part of the 
regularization procedure, we have introduced a cut off function. We found that in our model 
there exists a critical half-width, of the order of 1.3 GeV 2 , below which a non-perturbative 
solution can be obtained. Particularly noteworthy is that even though the dressed vertex 
near threshold is largely pseudovector in nature, we find sizable admixtures of pseudoscalar 
coupling at large off-shell nucleon momenta. 

It is important to realize that off-shell from factors and self-energies can not be directly 
measured. The observables in quantum field theory are obtained from the S-matrix, not 
the Green's functions. The latter will depend on the representation of the fields in the 
Lagrangian of the theory, whereas the former do not. In principle, observables should 
not depend on the regularization and renormalization procedures chosen. This, however, 
does not apply to the Green's functions. To draw qualitative conclusions about physical 
processes, one should treat all off-shell ingredients of a model consistently, calculating them 
with the same Lagrangian and adhering to the same model assumptions. For example, the 
present model for the pion-nucleon off-shell form factors and the nucleon self-energy can be 
consistently utilized in a K-matrix approach to ttN scattering, as shown in Section III C. 

Although the present work deals with nNN vertices where one of the nucleons is off-shell, 
we mention here that there exists a large amount of work on the form factor Gi(m 2 , m 2 , q 2 ) 
(see Eq. ([I])) for on-shell nucleons and an off-shell pion (see, e.g., |]32| , |35| , |36[| and references 



therein). In particular, in Ref. P2| a system of N, 7T, A, p, e, and u hadrons is considered 



in a consistent field theoretical framework. The approach of Ref. [35[] is based on a meson 



exchange model for ttN scattering, where, apart from the nucleon and the pion, also Delta 
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isobar and correlated 7nr exchange contributions are included. As mentioned before, the 
Delta resonance and the p meson are important ingredients in a quantitative description 
of ttN scattering and pion photoproduction. Therefore, one should expect these degrees of 
freedom to play a prominent role in a realistic model for pion-nucleon form factors. In the 
present work, only nucleon and pion fields are included, and only discontinuities associated 
with the one-pion threshold are taken into acoount. In this simplified model we focus on non- 
perturbative aspects of the consistent dressing of the pion-nucleon vertex and the nucleon 
self-energy. Contributions to the imaginary parts from higher thresholds can be included 
in our model either explicitely, by allowing intermediate states with two or more pions, 
or effectively, by considering baryon and meson resonances (for example, A and p). This 
extention of the model is in progress. 
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APPENDIX A: THE SELF-ENERGY 



Here some details on the evaluation of the imaginary part of the nucleon self-energy are 
given. To calculate the imaginary parts of the self-energy functions, we need to evaluate the 
pole contribution I po i e (p) as given in Eq. @. In general the integral can be expressed as 
Ipoieip) = l^ii.p) + h(p), where If and h are scalars in spinor space. Since the only Lorentz 
vector in the problem is If must be proportional to it, Ii(p) = ((h ■ p)/p 2 )p fl . Hence, 
one may write I po ie{p) — ^>h{p 2 ) + h(p 2 ), where h(p 2 ) = (h • p)/p 2 and hip 2 ) = h(p 2 ) are 
Lorentz scalars. They equal 



l x {p 2 ) = ? + 3 ^ 4 r(p 2 ) 6(p 2 - (m + p) 2 ) (Al) 



HP 2 ) = r(p 2 ) 6(p 2 - (m + p) 2 ). (A2) 



where r(p 2 ) = y \(p 2 , m 2 , p 2 ), with the Kallen function defined as\(x,y,z) = (x — y — z) 2 — 
Ayz. Using these expressions and introducing the shorthand notation, g 12 = ReGi 2 (p 2 ), the 
imaginary parts of the self-energy functions can be determined from Eqs. (|7],|§[3D, 

ImA^(p 2 ) = 3( - (<?! - g 2 ) 2 h( P 2 ) + 2( 9l - 92)9^- - gf^f 



m 77T 
2 1 ™2 ,,2 



9i 



5p 2 + m 2 -p 2 p 2 + m 2 -p 2 2 3p 2 + m 2 -p 2 x 
+ ( n + 2 )92 - 2 2 9i92 \ , (A3) 
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ImB + (p ) = 3 [{gi-g2) 2{g x -g 2 )g 2 — + g 2 — 

= 16^ r(p2)%2_(m + /i)2) 

„2 i o^2 _ ,,2 ^2 1 o^2 _ ,,2 



r 2 2p 2 + 2m 2 - /T 2 p 2 + 3m 2 - /i 2 1 . A ^ . 

x K + — 2 —92 - ~ 2 —9i92 \ , A4 



the factor 3 in the above equations results from the multiplication of the isospin matrices, 
T a T a = 3, and the minus sign in front of I\{p 2 ) from commuting the 7 5 matrices. The real 
parts of A n+l (p 2 ) and B n+l (p 2 ) are found by applying dispersion relations Eq. (TO), where 
all integrals are done numerically. 

For later use, the dressed nucleon propagator is written as 

(S(p))- 1 = Z 2 {j> - m B ) - [ReA(p 2 )i) + ReB(p 2 )m) = a(p 2 )(^ - £(p 2 )), (A5) 

where 

a( r?) = z 2 -R eA (r?), ^) = Z2(m - to) + feB(p2) '" . (A6) 

a(p z ) 

The spectral function of the self-energy is introduced through 

E L (p) = C(w)A + (fl + C(-w)A-(# , (A7) 

where A ± (p') = (±^ + u)/(2u) are the projectors on positive- and negative-energy states of 
the nucleon with the invariant mass u = \fp 1 ~ > 0. The spectral function can now be defined 

as 



T(±lu) = T-Im({±uj). (A8) 

7T 



Equating the right-hand side of Eq. ( |A7| ) and the form of £r,(p) from Eq. (|7|) yields 

T(±uj) = --(ujImA(p 2 )±mImB(p 2 )) . (A9) 

7T 



APPENDIX B: THE FORM FACTORS 

The calculation of the imaginary parts of the form factors can be reduced to computing 
one-dimensional integrals which are done numerically. First consider the integral on the 
right-hand side of Eq. (|12"D . J po i e can be split as J po i e = 7^1+ J2, where Jf and J 2 are scalars 
in spinor space, and a possible rank-2 tensor structure vanishes since u(p')fi = u(p')m. For 
the same reason Jf is proportional to only the vector p^. Following the same argumentation 
as used in Appendix A, we can write 

Jpoie = M + J 2 , (Bl) 

where we have introduced the Lorentz scalars J\ and J 2 . To write down expressions for J\ 
and J 2 , we define the following functionals: 
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Kilf] 



-i a{w' 2 Y 



-i 
i 



a(w 



r2\ 



dx ■ 



f(w' 2 ) 



-i «K)K-( 2 K))' 

-i "aK)K-( 2 K))' 



(B2) 
(B3) 
(B4) 
(B5) 



where / is any function for which the integrals exist, a and £ are given by Eq. (|A6|) for the 
(n+l) st iteration, and 



w 
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(jp' + k) = m 2 + fi 2 - 



2 p 4 — (m 2 — fj, 2 ) 2 



r(p 



2\2 



2p 2 



2p 2 



(B6) 



with x being the cosine of the polar angle between the three vectors ~p*' and k . Now J\ 
and J 2 can be written as 



J 1 = J 1 (p 2 ) = -U Kl -K, 



+(k 3 -k 4 



\p 2 + m 2 — fi 2 
V 



,p 2 + m 2 — fi 2 £ m 2 — 



Am 



2p 2 



i . £-4m 2 g\ 

{9192 + -^r 9 ^ + 2 



){mg 1 + (£- m)g 2 y 



(B7) 



and 



J2 = Mp 2 



+(k 3 + k. 

r(p 2 ^ 



+ 



, £ — 4m | p 2 + m 2 — fi 2 N ,, 



2m 

, £ — 2m p 2 + m 2 — /i 2 



+ 



Am 2 



)9i + £i#2 



2m 

% 2 -(m + /i) 2 ). 



Am 2 



) (m^i + (£ -m)g 2 y 



(B8) 



Since for a given / in Eqs. ( |B2| - p5|) the 1*Q are functions of p 2 only, Eqs. ( |B^ , |B8| , [BT|) show 
that Jpoie depends only on the Lorentz vector p and does not depend on p' as it might appear 
from the right-hand side of Eq. fll2|) . 



Using Eqs. (pl]Jil| ,|3|) , one obtains for the imaginary parts of the form factors: 
ImGTV) = 9iUP 2 ) + (^^2 - 9i) Ji (P 2 ) 



m z 



and 



ImG n 2 + \p 2 ) = g 2 J 2 (p 2 ) + (g 2 - g^p 1 



(B9) 



(BIO) 



where J\ and J 2 are given by Eqs. ( |B7| ) and (|B8|). Finally, Eq. fll3|) is applied to obtain the 
real parts of the (n+l) st iteration for the form factors. 



15 



REFERENCES 



[1] E. M. Nyman, Nucl. Phys. A160, 517 (1971). 

[2] S. Kondratyuk, G. Martinus, and 0. Scholten, Phys. Lett. B 418, 20 (1998). 
[3] Yi Li, M. K. Liou, and W. M. Schreiber, Phys. Rev. C 57, 507 (1998). 
[4] A. Yu. Korchin, O. Scholten, and F. de Jong, Phys. Lett. B 402, 1 (1997). 
[5] F. L. Gross, J. W. Van Orden, and K. Holinde, Phys. Rev. C 45, 2094 (1992). 
[6] F. Gross and Y. Surya, Phys. Rev. C 47, 703 (1993). 
[7] B. C. Pearce and B. K. Jennings, Nucl. Phys. A528, 655 (1991). 
[8] C. Schiitz, J. W. Durso, K. Holinde, and J. Speth, Phys. Rev. C 49, 2671 (1994). 
[9] T. Sato and T.-S. H. Lee, Phys. Rev. C 54, 2660 (1996). 
[10] K. Nakayama, A. Szczurek, C. Hanhart, J. Haidenbauer, and J. Speth, Phys. Rev. C 
57, 1580 (1998); J. Haidenbauer. C. Hanhart, J. Speth, K. Nakayama, and J. W. Durso, 
in Book of Abstracts of the 8-th International Conference on the Structure of Baryons 
"Baryons'98", Bonn, 1998, p. 84. 
[11] A. Bincer, Phys. Rev. 118, 855 (1960). 

[12] E. M. Nyman, Nucl. Phys. A154, 97 (1970); M. G. Hare, Ann. Phys. 74, 595 (1972). 

[13] H. W. L. Naus and J. H. Koch, Phys. Rev. C 36, 2459 (1987); P. C. Tiemeijer and J. 
A. Tjon, Phys. Rev. C 42, 599 (1990); J. W. Bos and J. H. Koch, Nucl. Phys. A563, 
539 (1993); H. C. Donges, M. Schafer, and U. Mosel, Phys. Rev. C 54, 950 (1995). 

[14] M. Ida, Phys. Rev. 136, B1767 (1964); G. N. Epstein, Phys. Lett. 79B, 195 (1978). 

[15] W. T. Nutt and C. M. Shakin, Phys. Rev. C 16, 1107 (1977); 

[16] N. N. Bogoliubov and D. V. Shirkov, Introduction to The Theory of Quantized Fields 

(Interscience Publishers, inc., New York, 1959), chapt. IX. 
[17] R. E. Cutkosky, J. Math. Phys. 1, 429 (1960); G. 't Hooft and M. J. G. Veltman, 

Diagrammar, CERN Yellow Report 73-09. 
[18] W. D. Brown, R. D. Puff, and L. Wilets, Phys. Rev. C 2, 331 (1970). 
[19] G. Krein, M. Nielsen, R. D. Puff, and L. Wilets, Phys. Rev. C 47, 2485 (1993); M. E. 

Bracco, A. Eiras, G. Krein, and L. Wilets, Phys. Rev. C 49, 1299 (1994). 
[20] C. A. da Rocha, G. Krein, and L. Wilets, Nucl. Phys. A616, 625 (1997). 
[21] P. F. A. Goudsmit, H. J. Leisi, E. Matsinos, B. L. Birbrair, and A. B. Gridnev, Nucl. 

Phys. A575, 673 (1994). 
[22] O. Scholten, A. Yu. Korchin, V. Pascalutsa, and D. Van Neck, Phys. Lett B 384, 13 

(1996). 

[23] S. Kamefuchi, L. O'Raifeartaigh, and A. Salam, Nucl. Phys. 28, 529 (1961). 

[24] S. Scherer and H. W. Fearing, Phys. Rev. C 51, 359 (1995); H. W. Fearing, Phys. Rev. 

Lett. 81, 758 (1998). 
[25] R. M. Davidson and G. I. Poulis, Phys. Rev. D 54, 2228 (1996). 
[26] J. Adam, Jr., F. Gross, and J. W. Van Orden, |nucl-th/9708"00l . 
[27] E. Kazes, Nuovo Cimento 13, 1226 (1959). 

[28] V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, Quantum Electrodynamics (Perg- 
amon Press, 1982), p. 474 and 476. 

S. Bauberger and M. Bohm, Nucl. Phys. B445, 25 (1995); A. Frink, J. G. Korner, and 
J. B. Tausk, |hep-ph/970949T?| ; P. Post and J. B. Tausk, Mod. Phys. Lett. A 11, 2115 
(1996). 



16 



[30] S. Weinberg, The Quantum Theory of Fields, vol. 1 and 2 (Cambridge University Press, 
1996). 

[31] A. Yu. Korchin, O. Scholten, and R. G. E. Timmermans, to be published in Phys. Lett. 
B. 

[32] J. Flender and M. F. Gari, Phys. Rev. C 51, R1619 (1995). 

[33] D. Drechsel, O. Hanstein, S. S. Kamalov, and L. Tiator, |nucl-th/9807001 . 

[34] Ulf-G. Meissner, Czech. Jour. Phys. 45, 153 (1995). 

[35] Schiitz, J. Haidenbauer, and K. Holinde, Phys. Rev. C 54, 1561 (1996). 

S. A. Coon and M. D. Scadron, Phys. Rev. C 23, 1150 (1981); A. W. Thomas and K. 
Holinde, Phys. Rev. Lett. 63, 2025 (1989); T.-S. H. Lee, |nucl-th/9502005| ; K.- F. Liu, 
S.-J. Dong, and T. Draper, Phys. Rev. Lett. 74, 2172 (1995); G. Holzwarth and R. 
Machleidt, Phys. Rev. C 55, 1088 (1997). 



17 



TABLES 



TABLE I. Values of the renormalization constants /, Z 2 and 5m for the two choices for the 
function G°(p 2 ), Eq. <M) and Eq. (M). 



Case 


/ 


z 2 


5m [MeV] 


(I) 


12.42 


0.848 


-57.4 


(II) 


12.43 


0.848 


-55.5 



18 



FIGURES 



FIG. 1. The graphical representation of the system of Eqs. (0). The thin and thick solid lines 
correspond to the free and dressed propagators of the nucleon, respectively. The dashed line is the 
propagator of the pion. The circle indicates the dressed ttNN vertex, and the square stands for the 
counterterm contribution to the nucleon self-energy. In the equation for the vertex the propagators 
of the external lines are stripped away, as indicated by the dashes on these lines. 



FIG. 2. The pole contribution of iteration n + 1 to the self-energy (the left picture) and the 
vertex (the right picture), as expressed by Eqs. (p[]l2]). The notation is as in Fig. 1, with the 
subscript R indicating that the vertex and propagator are calculated using only the real parts of 
the form factors and self-energy functions. The crosses on lines indicate that the corresponding 
particles are put on their mass shells. 
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FIG. 3. The imaginary parts of the pseudo-vector and pseudoscalar irNN form factors as func- 
tions of the momentum-squared of the off-shell nucleon, defined in Eq. (0). The calculations 
correspond to the two cut-off functions Eq. (|23|) (the left panels) and Eq. (|24]) (the right panels). 
The drawn (resp. dotted) curves show the converged (resp. first iteration) results. 
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FIG. 4. The real parts of the pseudovector and pseudoscalar ttNN form factors as functions 
of the momentum-squared of the off-shell nucleon, defined in Eq. (Q). The calculations correspond 
to the two cut-off functions Eq. ( ^3|) (the left panels) and Eq. (|||) (the right panels). In the top 
panels the zeroth iteration and the converged form factors are given by the dotted and drawn lines, 
respectively. In the middle panels the converged results and those of the first iteration are shown 
for the pseudoscalar form factor and the loop contribution to the pseudovector form factor. The 
bottom panels show the ratio of the pseudoscalar and pseudovector form factors, where the drawn 
(resp. dashed) curves correspond to the converged (resp. first iteration) results. 



FIG. 5. The imaginary parts of the self-energy functions A{p 2 ) and B(p 2 ), as defined in Eq. ([?]). 
The calculations correspond to the two cut-off functions Eq. ( |23|) (the left panels) and Eq. ( p4|) (the 
right panels). The drawn (resp. dotted) curves are the converged (resp. first iteration) results. 
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FIG. 6. The same as in Fig. 5, but for the real parts of A(p 2 ) and B(p 2 ). 
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FIG. 7. The self-energy spectral function T(uS) as function of invariant mass of the nucleon. 
The upper (resp. lower) panel corresponds to the calculations with the cut-off function Eq. (23) 
(resp. Eq. (|2"4|)). The drawn (resp. dashed) curves are the converged (resp. first iteration) results. 
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